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Abstract 

We reconsider the phenomenological implications of two texture zeros in sym- 
metric neutrino mass matrices in the light of the recent T2K results for the reactor 
angle and the new global analysis which gives also best fit values for the Dirac CP 
phase 5. The most important results of the analysis are: Among the viable cases 
classified by Frampton et al. only Ai and A2 predict #13 to be different from zero 
at 3a. Furthermore these two cases are compatible only with a normal mass spec- 
trum in the allowed region for the reactor angle. At the best fit value Ai and A2 
predict 0.024 > sin 2 13 > 0.012 and 0.014 < sin 2 13 < 0.032, respectively, where 
the bounds on the right and the left correspond to cos 5 = — 1 and cos 5 = 1, re- 
spectively. The cases Bx, B2, B3 and B4 predict nearly maximal CP violation, i.e. 
cos <5«0. 

PACS-numbers: 14.60.-z, 14.60.Pq, 14.60. St, 23.40.Bw. 

1 Introduction 

Recently the T2K Collaboration pQ gave hints for a nonzero reactor angle, and also the 
results of the MINOS Collaboration |2J point towards the same direction. The global fits 
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of neutrino oscillation experiments giv^ 

0.001 < sin 2 # 13 < 0.035 (NS), 

0.00K sin 2 9 13 < 0.039 (IS), U 

K h (1) 

0.005 < sin 2 9 13 < 0.050, $\ 

and the best fit values are sin 2 #13 = 0.013 and sin 2 #13 = 0.025, respectively. In particular 
the global analysis by Schwetz, Tortola and Valle [3] gives a weak hint for a nonvanishing 
CP violating phase, namely (at the best fit point) 

sin 2 13 = 0.013, <5 = -0.61tt (NS), (2) 
sin 2 #i3 = 0.016, 5 = -0.41tt (IS). (3) 

A lot of papers have been proposed recently in order to reproduce such a large value 
of the reactor mixing angle [5]. Already before the recent T2K data there have been 
models based on discrete flavor symmetries which predict a large reactor mixing angle — 
for an incomplete list see Ref. [6], and for a classification of models with flavor symmetries 
classified by their predictions for the reactor angle see [7]. 

Here we reconsider the interesting case of Majorana neutrino mass matrices with two 
texture zeros in the basis where the charged lepton mass matrix is diagonal, which has 
been extensively studied in the past. Our aim is to point out the phenomenological 
implication of such textures in the light of the T2K result^. 

It was shown in [TQl [11] and [8] that, in the basis where the charged lepton mass 
matrix is diagonal, there are seven types of two texture zeros in symmetric neutrino mass 
matrices compatible with the experimental data on neutrino oscillations. In this work we 
want to analyze the correlation between the CP violating phase 5 and the reactor mixing 
angle #13 in the framework of these two texture zeros. 

Another interesting possibility is to place texture zeros in the inverted neutrino mass 
matrix-see e.g. [12J. The implications of this type of two texture zeros on the reactor 
mixing angle and CP violation have been studied in [13J. 

2 Two texture zeros 

In the basis where the charged lepton mass matrix is diagonal, the lepton mixing matrix 
U and the Majorana neutrino mass matrix Ai u are related via 

M v = U* diag(mi, m 2 , m 3 ) U ] . (4) 

The standard parameterization [H] for U is given bj|f| 

U = e i& Ve ia , (5) 



1 Throughout this work the abbreviations NS and IS will stand for normal and inverted neutrino mass 
spectrum, respectively. 

2 While we were finishing this work two papers treating the same problem were published in [8j|9]. 
3 The parameterization used here is a re- writing of the symmetrical parameterization proposed in |15j . 
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where 



/ 



V 



C13C12 



C13S12 



s 13 e- iS \ 



-C23S12 - S23Sl3Cl 2 e 
AS 



1 5 



C23C12 - S 2 3Sl3Si 2 e' u S23C13 
iS 



(6) 



\ S23S12 — C23Si3Ci2e to — S23C12 — C23Si3Si2e C23C13 J 

a = diag(ai, a 2 , a 3 ) and a = diag(<7i, cr 2 , cr 3 ). Cij = cos^-, = sinfly with 8^ G [0, 7r/2]. 
<5 G [0, 2tt) is the CP violating phase and the <Tj G [0, 2tt) are the Majorana phases, which 
are not measurable in oscillation experiments. The phases oij are irrelevant for neutrino 
oscillations and will play no role in our analysis, as we will see in the following. Inserting 
([5]) into ([1J we obtain 

( M "Yij = J2 mkU ^ k = J2 mke2mkV ^ kel{at+a]) - ( ? ) 

k k 

Placing a texture zero in the neutrino mass matrix corresponds to the condition 

{M v )ii = U (^(^ = 0) (8) 
for some indices (i, j). Defining := mke 2%<Jk and dividing by e l ^ ai+a ^ we arrive at 

^2fikVikV jk = 0. (9) 
k 

The assumption of two texture zeros can thus be described by the two equations 

y~l HkVakVbk = 0) /] HkVckVdk = 0- (10) 



The viable cases of two texture zeros given in pU] and the corresponding parameters 
(a, 6, c, d) can be found in table [TJ 



case 


texture zeros 


(a,b,c,d) 


Ai 


(M v ) ee = (M u ) efl 


= 


(1,1,1,2) 


A 2 


(M v ) ee = {M u ) er 


= 


(1,1,1,3) 


Bi 


(M v )w = {M v )er 


= 


(2,2,1,3) 


B 2 


{M V ) TT = (M v )ep 


= 


(3,3,1,2) 


B 3 


(M v )w = {M v ) etl 


= 


(2,2,1,2) 


B 4 


(M v )tt = (M u ) er 


= 


(3,3,1,3) 


C 


{M v )w, = {M V ) TT 


= 


(2,2,3,3) 



Table 1: The viable cases in the framework of two texture zeros in the Majorana neutrino 
mass matrix Ai v and a diagonal charged-lepton mass matrix Aii |10j . 
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3 General remarks 



The system (TlOT) is equivalent to 



14lVfcl V a 2Vb2 
VdVdl V c 2Vd2 




"M3 



V a3 V b3 



y c3 v d3/ 

The set of solutions of this system of linear equations depends on the determinant 



D 



abed 



det 



For D abcd ^ we find 




V cl V dl 



M3 



K 2 ^2 



V a lVblV c 2Vd2 — V a 2Vb2V c iVdl. 



(12) 



abed 



Vc2Vd2 

-VciVdi 



-v a2 v b2 
v al v bl 



v a3 v b3 
v c3 v d3 



(13) 



Since at least two neutrino masses must be nonzero, the above equation implies that the 
lightest neutrino mass is different from zerofl Thus we are allowed to divide by jj 3 and 
we can easily calculate 



r :- 





2 




2 














A»3 





(14) 



Inserting ( I13p into (j!4p we find an equation which relates the six quantities 

^31) 6*12, #23, dl3 



Ato2!, AmL 0i2, 6*2,s, 6*i3 and 5. 



Fixing the mass squared differences and the two mixing angles 9u and 623 (e.g. to their 
best fit values or their ncr-ranges), we obtain a relation between the reactor mixing angle 
#13 and 5. Note that in this way one can eliminate the unknown absolute neutrino mass 
scale. This approach has been previously used in [16J. 

The main question we have to answer before beginning our analysis is whether the 
determinant D abcd can become zero for the seven different cases within the experimental 
limits. The first issue we notice is that all entries of the 2 x 2— matrix 



KiVbi K2H2 

v cl v dl V c2 V d2 



(15) 



are nonzero (by experiment). Thus D abcd = 

ValVbl 

v cl v c 



dl 



implies 

V a 2V b 2 
V c 2V d 2 



(16) 



In general a normal (inverted) neutrino mass spectrum allows /ii = = 0). However, one can 
verify that within the experimental 3<r-range (|13j) implies fix = ^ [13 = for all types of two texture 
zeros we will study in this work. Thus the lightest neutrino mass must be nonzero. 
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Using the experimentally known fact that the absolute values of all elements of the second 
column of V are of the same order of magnitude, we find 







v a2 v b2 


v cl v dl 




V c 2Vd2 



\v al v t 



bl\ 



\v cl v c 



(II I 



(17) 



From 



\V a \ ~ (2/VE, 1/V6, l/v / 6) r 



'It 



one easily finds that the only case allowing ( ITT)) is C. Therefore for the cases Aj to B 4 we 
can assume -D a fecd 7^ and use equation (fT4"I) for our analysis. 

Let us now turn to case C. In [18] it has been shown that for 8 is = the determinant 
-D2233 becomes zero and the system ( TTOT) is therefore singular in this case. Inversely 
assuming -D2233 = 0, we can proceed as follows. Defining e = si^e 15 we find 



D 



2233 



sin(2# 12 )sin(2fl 23 )(l 



- ecos(20i 2 )cos(20 : 



23, 



Thus -D2233 can be zero only for e = or 

tan 2#i2 tan 2^23 



2e 



1 + e 2 



For < s? 3 < 0.05 we find 



2e 



1 + e 2 



2s 13 



1 1 + sf 3 exp(2z<5)| 



< 



2s 



13 



1 



< 0.48. 



'13 



Using the 3a— ranges provided in [3] one easily finds that at 3cr 

tan 20i2 tan 2023 > 8.55, 



(19) 



(20) 



(21) 



(22) 



which implies that S13 = is indeed the only possibility for -D2233 to become at 3cr. 
Since we are not interested in the limit S13 — > 0, we can use (ITU)) and (IT41 also to analyze 
case C. 



Analysis of the relation between sin 2 13 and cos 5 

It turns out that for all texture zeros studied in this work r (see equation (EH))) can be 
expressed as a rational function of at most cubic polynomials in cos 5, i.e. 

r = « (23) 
g(cos 0) 

where p and g are polynomials of order at most 3. Thus we find 

r g(cos<5) — p(cos5) = 0, (24) 

which is an equation of at most third order in cos 5. Thus the dependence of cos 5 on the 
mixing angles can be computed exactly. Note that (124")) will in general have more solutions 
than ( l23l) . because we have multiplied by q(cosS). In fact we have the additional solution 

g(coscT) =p(cos5) = 0, (25) 
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which corresponds to the limit Am^/m 2 , — > (see equation (TT4"j) ). i.e. a quasi-degenerate 
neutrino mass spectrum. 

For the cases Ai and A 2 fl2l]) is linear in cos 5. B 1; B 2 and C lead to quadratic equations 
and B 3 and B 4 yield cubic equations for cos 5, respectively. We used Mathematica to 
obtain the coefficients of ( |24|) . The well-known formulae for the general solutions of 
quadratic and cubic equations were implemented in C-programs which, scanning over 
the experimentally allowed ranges for r, # 12 and # 23 , allowed us to plot sin 2 # 13 versus 
cos S. sin 2 # 13 was varied between and 0.05 and for all other experimentally accessible 
quantities we used the values obtained from the newest global fit [3j including already the 
new T2K data pp. Our numerical analysis consists of the following steps: 

— Input: sin 2 # 13 , sin# 12 , sin# 23 , r = Am^/Amjj (best fit values or no— range (n = 
1,2,3)). The range of sin 2 13 (0 to 0.05) is divided into 600 steps. The ranges of 
r, sin 9 2 3 and sin # 12 are divided into steps of equal length in such a way that the 
corresponding la— ranges are divided into 40 steps. Thus e.g. the computation for 
the lcr-range alone consists of 600 x 40 3 = 3.84 x 10 7 individual calculation cycles. 

— Equation (I2~4j) is solved for cos 5 (there can be up to three solutions). 

— Only the real solutions G [—1,1] are processed further, the others are discarded. 

— Now we want to insert the remaining solutions for cos 5 into ( TT3|) to calculate the 
mass ratios 

m _ \m_ 

In order to do that we have to calculate e iS from cos 5. There are two solutions to 
this problem, namely 

e iS ■= cos5±i y/1 - cos 2 5, (27) 



(26) 



but since they are complex conjugates of each other the mass ratios (1261) do not 
depend on the choice of the solution. 

— Finally the program checks whether the following inequalities are fulfilled. 

NS), 

(28) 



mi 


< 1, 


m 2 


<1, 


mi 


< 1 


m 3 




m 3 




m 2 




mi 


>1, 


m 2 


>1, 


mi 


< 1 


m 3 




m 3 




m 2 





If they are fulfilled the data point (cos 5, sin 2 13 ) is stored. 



— Finally, when the whole parameter range has been scanned, all stored data points 
are plotted^- 



5 In order to create plots of a suitable size (in terms of disk space) we constructed a lattice dividing 
the range of sin 2 #13 into 600 and the range of cos<5 (—1 to 1) into 800 points. Data points falling into 
the same part of the lattice were plotted only once. 
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It turns out that a simple scan over the allowed ncr-ranges for the parameters as 
described above works very well for all cases of types A and B. However, for case C (NS) 
the sizes of the steps chosen in our systematic scan are just too large in order to obtain 
enough data points to produce good and reliable results. The reason for this issue is that 
C (NS) implies #23 ~ 45° [IB] , so one would need an enormously high resolution in the 
scan over the ncr-ranges of sin #23 to produce reliable results. Thus we have to analyze C 
(NS) in a different way. A good method to deal with case C (NS) is to assign the input 
parameters sin 2 # 13 , sin# 12 , sin #23, r = Am^/Am^ random values in their ncr-ranges, 
rather than varying them step by step. In this way one obtains a so-called scatter plot. 
Since also case C (IS) shows some hints of problems using a systematic scan, we also did 
a scatter plot for this case. The number of random points (sin 2 #13, sin 6*12, sin #23, r) we 
used was 10 9 for each of the ncr-ranges (n = 1, 2, 3). 

4 Results 

We will now present the results of our numerical analysis. As already explained we have 
produced plots of sin 2 9 13 versus cos 6 (see figures [H 4T2]) . The color code is the same for 
all plots: 

• The best fit value for the point (cos 6, sin 2 6>i 3 ) according to the global fit [3] is 
indicated by a black cross. 

• The best fit values for sin 2 #13 as a function of cos 6 according to our analysis are 
indicated by a black line. 

• The na— regions are shown as colored areas in the plots (red=lcr, green=2cr and 
blue=3cr). 

The cases Ai and A2 are incompatible with an inverted spectrum if the reactor angle is 
varied within the 3cr-range < s 2 3 < 0.05. Assuming a normal neutrino mass spectrum 
Ai and A 2 predict #13 to be different from zero at 3cr. For the best fit values for the 
observables, namely 9 12 , 2 3 and r, we predict 0.024 > sin 2 (9 13 > 0.012 corresponding to 
the bounds —1 < cose; < 1 for the case Ai and 0.014 < sin 2 #13 < 0.032 corresponding to 
— 1 < cos 5 < 1 for the case A 2 . 

The cases Bi, B2, B3 and B4 predict the Dirac CP phase to be close to maximal, i.e. 
cos 5 ~ 0. Note furthermore that the cases Bi (NS), B 2 (IS), B 3 (NS) and B 4 (IS) are 
incompatible with #23 > 45° .0 Therefore, for these cases, the plots do not show the black 
"best fit line" (the best fit value for sin 2 #23 given in [3] is 0.52 which corresponds to an 
atmospheric mixing angle larger than 45°). The generic predictions we found for the cases 
of type B concerning the atmospheric angle are shown in table [2J 

For case C and an inverted spectrum we do not find a strong correlation between 
sin 2 9 ri and cos S. However, from figure [12] we can see that the best fit value for the point 
(cos S, sin 2 612) lies in the lcr-region of the plot. For case C and a normal spectrum there 
is no correlation between the Dirac CP phase and the reactor angle (see figure fTTi) . but, 

6 This is in accordance with the best fit results for sin 2 #23 given in [17] for B3 and B4. 
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case 


NS 


IS 


Bi 


e 23 < 45° 


023 


> 45° 


B 2 


023 > 45° 


023 


< 45° 


B 3 


023 < 45° 


023 


> 45° 


B 4 


023 > 45° 


023 


< 45° 



Table 2: Inequalities for the atmospheric mixing angle for the cases of type B. 



CD 
C 




-1 -0.5 0.5 1 

cos 8 

Figure 1: The relation between sin 2 #13 and cos S for case Ai (normal spectrum). For 
description of the colors see the text. 



as was pointed out by Grimus and Lavoura [18], atmospheric neutrino mixing is close to 
maximal. As shown in figure H~3l there is a correlation between the reactor angle and the 
atmospheric angle, but the deviation from the maximal value of the atmospheric angle is 
negligible. 



5 Conclusions 

In the light of the recent T2K result which points towards a large reactor mixing angle 
#13, we reconsidered the interesting case of two texture zeros in the neutrino mass matrix. 
In particular we studied the correlation between the reactor mixing angle #13 and the 
Dirac CP phase S for the viable cases classified in (TO] as Ai, A2, Bi, B2, B3, B4 and C. 
All of these cases are still compatible with the global fit of the neutrino data at 3a, but 
only the cases Ai and A 2 predict the reactor angle to be different from zero at 3cr. In 



8 




Figure 2: The relation between sin 2 # 13 and cos 5 for case A 2 (normal spectrum) 



<x> 



0.05 



0.04 



0.03 



0.02 



0.01 




-1 



-0.5 



0.5 



cos 8 



Figure 3: The relation between sin 2 # 13 and cos 5 for case B x (normal spectrum) 
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Figure 4: The relation between sin 2 # 13 and cos 5 for case Bx (inverted spectrum) 
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Figure 5: The relation between sin 2 # 13 and cos 5 for case B 2 (normal spectrum) 
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Figure 6: The relation between sin 2 # 13 and cos 5 for case B 2 (inverted spectrum) 
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Figure 7: The relation between sin 2 # 13 and cos 5 for case B 3 (normal spectrum) 
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Figure 8: The relation between sin 2 # 13 and cos 5 for case B 3 (inverted spectrum) 



CD 

c 



0.05 



0.04 



0.03 



0.02 



0.01 




Figure 9: The relation between sin 2 # 13 and cos 5 for case B 4 (normal spectrum) 
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ure 10: The relation between sin 2 # 13 and cos<5 for case B 4 (inverted spectrum 
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Figure 12: The relation between sin 2 # 13 and cos<5 for case C (inverted spectrum) 
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gure 13: The relation between sin 2 6 , 13 and sin 2 #23 for case C (NS) (scatter plot for the 
-range with 10 7 random points). 
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particular for the case A 1; asserting all the free parameters their best fit values, predicts 
0.012 < sin 2 (9 13 < 0.024 while for the case A 2 assuming the best fit values predicts 
0.014 < sin 2 ^, < 0.032. 
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